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Loop quantum cosmology is shown to provide both the dynamical law and initial conditions for 
the wave function of a universe by one discrete evolution equation. Accompanied by the condition 
that semiclassical behavior is obtained at large volume, a unique wave function is predicted. 
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Traditionally, physical systems are modeled mathe- 
matically by providing laws governing the dynamical be- 
havior and specifying initial (or boundary) conditions. 
The latter select a particular solution to the laws, but 
usually all of them are allowed and describe the sys- 
tem under different conditions. However, in cosmology 
the situation is different: there is only one universe, and 
therefore only one fixed set of initial conditions can lead 
to the physically realized situation, fn this context, the 
big bang singularity is regarded as the point of "creation" 
of the universe at which initial conditions (or equivalent 
restricting requirements) have to be imposed. Since grav- 
ity is strong at that stage and classical general relativity 
breaks down (signaled by the appearance of a classical 
singularity), a quantization of the gravitational field is 
needed bringing us in the realm of quantum cosmology. 

The standard approach to quantum cosmology consists 
in quantizing a minisuperspace model which is obtained 
by specifying symmetry conditions, usually homogene- 
ity and isotropy, for the allowed metrics in space-like 
slices of a universe. This reduces the infinitely many 
degrees of freedom of general relativity to finitely many 
ones allowing standard quantum mechanical methods 
jl],^| . Due to general covariance the dynamical law is pro- 
vided by a constraint equation which takes the form of a 
second order differential equation — the Wheeler-DeWitt 
equation — for the wave function ip(a,4>) depending on the 
scale factor a > (conventionally used as internal time) 
and matter degrees of freedom collectively denoted by <f>. 
However, the classical singularity remains, and no initial 
conditions are provided by the formalism which leads at 
least to a two-parameter family (not counting matter de- 
grees of freedom) of solutions and not a unique (up to 
norm) one. The original hopes [jjj that there might be a 
unique solution to the constraint equation are not real- 
ized. 

To address this issue, proposals have been developed 
by several authors. However, these proposals have con- 
siderable arbitrariness since they are driven primarily by 
the authors' intuition as to how the classical singularity 
might be smoothed out by quantum gravity. Most well- 
known are the "no-boundary" proposal ||] and the "tun- 
neling" proposal which both describe the "creation" 
of a universe at the place of the classical singularity. In 
all those approaches matter is regarded as being irrele- 



vant in the early stages, and so the wave function ip is 
assumed (implicitly or explicitly ||) to be independent 
of the variables <f> for small a; this is already an initial 
condition which strongly restricts the ^-dependence of ip. 
We will take the same point of view concerning matter 
degrees of freedom here, which we regard as being justi- 
fied thanks to the dominance of gravity in early stages of 
the evolution. 

But still, there is a two-parameter family of solutions 
yj(a) from which one parameter has to be fixed (since the 
norm is irrelevant). This not only influences the wave 
function close to the singularity, but also its late time 
behavior because it selects a particular linear combina- 
tion of the expanding and contracting components in a 
WKB-approximation. However, as an initial condition 
it is specified at the classical singularity (e.g., by fixing 
the value "0(0) or by introducing an ad hoc "Planck 
potential" ||), and thus involves Planck scale physics for 
which we need a full quantum theory of gravity. 

One candidate for a quantization of general relativity 
is quantum geometry (see e.g. 0J|) which predicts dis- 
crete eigenvalues of geometrical operators like area and 
volume [9|-pl|. A symmetry reduction [[l2| of the quan- 
tized (kinematical) theory to cosmological models leads 
to loop quantum cosmology [[L3[, in which the discrete- 



ness of the volume is preserved ]f4| . All techniques used 
in this framework of quantum cosmology are very close 
to those of full quantum gravity, in contrast to standard 
quantum cosmology which is based on a classical symme- 
try reduction to a simple mechanical system and subse- 
quent quantization. Hence, the results of loop quantum 
cosmology should be more reliable, in particular close to 
the classical singularity where the two approaches show 
the largest differences. In fact, loop quantum cosmol- 
ogy has a discrete evolution equation |l5| , [l6|] which re- 
places the Wheeler-DeWitt equation and is singularity- 
free fr^-]f9||. The fact that the Hamiltonian constraint 
operator of loop quantum cosmology ff5|| is very close to 
that of the full theory |^(| gives rise to the hope that the 
results of JlTj can be extended to less symmetric models. 

Loop and standard quantum cosmology deviate most 
when applied right at the classical singularity. In this let- 
ter we will show that the particular form of the evolution 
equation of loop quantum cosmology, applied at vanish- 
ing scale factor, leads to a consistency condition for the 
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initial data. In this way the evolution equation provides 
both the dynamical law and initial conditions: dynamics 
dictates the initial conditions. Accompanied by a classi- 
cality condition for the solutions, a unique (up to norm) 
wave function is predicted. 

Isotropic loop quantum cosmology. In the triad rep- 
resentation of isotropic loop quantum cosmology [p~S| ] 
the scale factor a e M+ is replaced by a discrete label 
n € Z which parameterizes eigenvalues of the triad op- 
erator. An orthonormal basis of the kinematical Hilbcrt 
space is given by quantum states \n) labeled by the triad 
eigenvalue n which also determines volume eigenvalues: 
y\ n ) = V(\n\-X)/vt\n) with 

^• = (7i|)V^' + 5)0' + l) (1) 

(7 G M + , which is of order one, is the Barbero-Immirzi 
parameter labeling inequivalent representations of the 
classical Poisson algebra, and lp — \fnh with re = 8ttG is 
the Planck length). The volume operator has eigenvalue 
zero with threefold degeneracy (for the states | ± 1) and 
|0)), but only one of them, |0), has degenerate triad and 
so corresponds to the classically singular state. The wave 
function ijj(a,4>) of standard quantum cosmology is re- 
placed by the coefficients s„(</>) of a state |s) = Y^, n s n\n). 
For large \n\ the correspondence between a and n is 
\n(a)\ ~ 6a 2 7 _1 /p 2 which follows from the volume spec- 
trum (\n\ = 2j + 1). 

The Hamiltonian constraint equation for spatially flat 
models takes the form of a discrete evolution equation 
(see [[l9| for the case of models with positive spatial cur- 
vature): 

i(l + 7 - 2 )4 8) Sn+8W>) - A^s n+i {cj>) - 24°) S „(0) 

-4r 4) ^- 4 (</>) + \{i + 7 - 2 )^r 8) ^-8^) 

= -\ 1 KllH^{n)s n {4>) (2) 

with 

4 ±8) : = (Vjn±8|/ 2 ~ Vj n±8 | /a _0 *£fc8*»±4 (3) 
4 ±4) := (V|n±4|/2 - Vj n±4 | /a _i) (4) 
4 0) == (V|n|/2 " V\n\/2-l) 

x (|(1 + !- 2 ){k-k+ +i + k+k-_ 4 ) - 1) (5) 

where Vj are the eigenvalues (Q) of the volume operator 
with V-i — 0, and the coefficients can be chosen to be 
non- vanishing by a suitable ordering of the extrinsic cur- 
vature operator and are approximately sgn(rt) for large 
\n\ (see [ p^[ for explicit expressions in terms of the volume 
eigenvalues). Here we introduced a matter Hamiltonian 
whose particular form is irrelevant. It only matters 
that it acts diagonally in the triad degrees of freedom 
which is always the case in the absence of curvature cou- 
plings. 



Pre-Classicality. Compared to the standard 
Wheeler-DeWitt equation of second order the discrete 
evolution equation is of order sixteen. So the problem 
of a unique solution seems to be more severe at first 
sight, but many of the additional solutions can easily be 
seen to not correspond to a semiclassical solution. Let 
us call a wave function s n pre-classical if and only if, 
at large volume (n 3> 1), it is not strongly varying at 
the Planck scale (increasing the large label n by one), 
although it may oscillate on much larger scales (increas- 
ing n by a macroscopic amount). Note that a difference 
equation with fixed step size, as is always the case here, 
may have solutions which are very different from those 
of an approximating differential equation even though 
all solutions of the differential equation are well approx- 
imated when the step size goes to zero ^J. A common 
possibility is a solution with alternating sign between 
successive n, which cannot correspond to a continuous 
solution of a differential equation. Due to instabilities 
there can also be solutions with exponentially increasing 
absolute value, even in regimes where the solutions of 
the differential equation are purely oscillating (i.e. in the 
classically allowed range in a WKB approximation). A 
precise formulation of the phrase "not strongly varying" 
can be given in the following way. Note first that the 
Barbero-Immirzi parameter 7 enters a(n) — ~flpn/6, 
which is used here as internal time. Although the phys- 
ical value of 7 is fixed and of order one [plj , we can use 
the 7 — > limit, together with n —> 00 such that a(n) is 
fixed, to decide whether a wave function is pre-classical. 
In this limit the difference a(n + 1) — a(n) becomes in- 
finitesimal implying a continuum limit. A wave function 
s n is pre-classical if and only if its limit 7 — > 0, n — > 00 
exists, providing a rigorous check of the pre-classicality 
condition. Note that re and h, and so lp , are fixed in this 
limit and we are still dealing with quantum cosmology. 
In fact, standard quantum cosmology can be shown to 
be the above limit of loop quantum cosmology. 

Our condition picks out only those solutions which are 
oscillatory on large scales but almost constant on the 
Planck scale. Since this is a pre-requisite for a sub- 
sequent WKB-approximation, we call it pre-classicality. 
Whenever it is fulfilled, a discrete wave function s n ((f>) 
can be approximated at large n by a standard con- 
tinuous wave function tp(a) :— s n / a ) with n(a) = 
6a 2 7 _1 /p 2 as above, which approximately solves the 
standard Wheeler-DeWitt equation up to corrections of 
order y/jlp/a p7| , |l9| |. Thus, standard quantum cosmol- 
ogy is realized only as an approximation valid at large 
volume where the discreteness of quantum geometry is 
irrelevant (see Fig. 0). 

Since the Wheeler-DeWitt equation is of second order 
and so has two independent solutions, there can be at 
most two independent pre-classical solutions s of our 
discrete evolution equations, such that any pre-classical 
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solution can be written as s = as + + bs with a, b G C. 
To demonstrate this explicitly, we introduce 

tm '■= 7 _1 ^p 2 (^2|m| - ^2|m|-l) s 4m 

P(m) := i7K^i?0(m)(V r 2 | m | - V 2 \ m \-iy 1 , 

using the expectation value H^fn) of H^n) in a matter 
state, such that for \n\ ^S> 1, where fc„ ~ sgn(n), the 
evolution equation (|2|) takes the form 

|(1 + j- 2 )t m+2 - t rn+1 + (|(3 - 7 ~ 2 ) + P{m))t m 
-t m _i + i(l+ 7 - 2 )< m _2 = 0. (6) 

In a classical regime \m\ is large and P(m) ~ ^nH^/a 
is approximately constant on a range small compared to 
|m|. In this case we have a linear difference equation 
with constant coefficients whose solutions can be found 
by an ansatz t m cx e m0 with € C which in (||) yields 
the quadratic equation 

(1 + 7~ 2 ) cos 2 - 2 cos + 1 - 7~ 2 + P = 

which has solutions 

cos0 = (1 + 7" 2 )- 1 (l± -{l+l- 2 )P) 

being real with modulus smaller than one such that is 
real when 7 is of the order one and P is small. 

If the matter does not contribute a Planck size en- 
ergy, P is small and we have cos0o = 1 — e + 0(e 2 ) 
or cos^i = (1 + 7~ 2 )~ 1 ( 1 - 7 ~ 2 ) + e + 0{e 2 ) with 

< e := \^ 2 P <C 1. The first possibility, expand- 
ing cos $o = 1 — 5^0 + 0(0 0)1 leads to two solutions 
&o = jVP~+0(P) and -0 O with |0 O | < 1 both of which 
imply pre-classical t m = e ±lme ° . Because 7 is not large 
compared to one, the second possibility cos0i leads to 

01 which violates pre-classicality (e.g., for 7 = 1 we have 
01 = ±7r/2 and t m oc (±i) m ). 

All 16 independent solutions t n /i = "/~ 1 lp' 2 (V\ n \/2 — 
V\ n \/2-i)s n of (||) can be obtained as 

. _ ±™0 o /4 ±in0i /4 
! -n/4 — e 1 e > 

/_2^n e ±in0o/4 (_2)7J e ±in.0i/4 
cr n e ±in0 o /4 Qr o .n e ±in0i/4 

where cr can be +« or — i. Obviously, only the first two 
are pre-classical. The definition using 7 — > (o finite) is 
applied as follows: with n = 6a 2 7 _1 Zp 2 we have 

lini 7" 1 Zp 2 (V|„| / 2 - VLi/2_i) = a/2 , 

lim e ±inf> °/ 4 = exp (±§iVPa 2 /Z 2 ) , 

whereas lim 7 ^o 0i = tt and so the limit for 7 — > 0, n — > 00 
does not exist for the remaining 14 solutions. 



Thus, by using the pre-classicality requirement, which 
is a prerequisite for any semiclassical analysis, we arrive 
at the same situation as in standard quantum cosmology: 
there are two independent solutions from which we have 
to select a linear combination up to norm. 
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FIG. 1. The unique solution (+) of the Hamiltonian con- 
straint (Q), when the Ricci curvature is only due to a pos- 
itive cosmological constant A := IpA = 2 • 10~ 4 (7 = 1), 
which is pre-classical for large positive n. Evolving backwards 
through n = 0, the wave function picks up a wildly oscillating 
component and is no longer exactly pre-classical at negative 
n. The standard quantum cosmology wave function tp(a), 
subject to (i£(4a) _1 d/da(a _1 d/da) + 3Aa 2 Z~ 2 )y^/>(a) = 
in the ordering corresponding to (^|), is given by 
ip(a) = a'i(AAi(~(3X)ia 2 lp 2 ) + BBi(~(3X)ia 2 lp 2 )) in 
terms of Airy functions. Wave functions i/j(a) for two choices 
of the parameters A and B are shown: the continuous line is 
the unique (up to norm) choice fulfilling DeWitt's V>(0) = 0, 
which in this case is in good agreement with s n at positive n, 
whereas any other choice leads to a diverging wave function 
(dashed line). 

Dynamical Initial Conditions. Up to now we consid- 
ered only the semiclassical regime, but there is an addi- 
tional feature of loop quantum cosmology which emerges 
right at the classical singularity, deeply in the Planck 
regime where the approximation by standard quantum 
cosmology breaks down: the highest order (or lowest or- 
der when we evolve backwards) coefficient vanishes when 
we try to determine sq. At first sight, it seems that this 
is a breakdown of the evolution similar to the classical 
situation. However, as demonstrated in ]l7|,[lj|, this is 
not the case in the particular factor ordering of the con- 
straint chosen above because so completely drops out of 
the evolution equation. (This observation depends cru- 
cially on the fact that ff^so(</>) = which is always true 
in quantum geometry p^Jl9[ | but would be impossible 
without space-time discreteness.) 

Instead of determining so the evolution equation leads 
to a consistency condition for the initial data: starting 
from a general pre-classical solution s n = as„ + bs~ for 
large n and evolving backwards, we eventually arrive at a 
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point where we have to apply for n — 8. At this value 

( — 8) 

of n the lowest order coefficient A y s ' vanishes as noted 
above causing s to drop out. Since by assumption we 
have already determined all s n for n > (for which there 
is no vanishing coefficient in the evolution equation), the 
would-be equation for sq leads to a further condition for 
higher s n (34, sg, S12 and s±e) which upon inserting the 
general pre-classical solution s n = as„ + bs~ implicitly 
yields a linear relation between the two free parameters 
a and b. This leaves us with a unique solution (up to 
norm) . 

Conclusions. We have shown that loop quantum 
cosmology implies a discrete evolution equation which 
uniquely determines a state (up to norm) behaving semi- 
classically at large volume. It is important to adapt the 
standard condition for semiclassicality in a WKB approx- 
imation taking the discreteness of time into account. This 
leads already to a strong reduction of the allowed solu- 
tions, but the crucial condition for the uniqueness arises 
only from the particular structure of the evolution equa- 
tion in quantum geometry. We remark that in general 
it is only possible to require pre-classicality at one con- 
nected domain of large volume. If one evolves through 
a classical singularity, the wave function may pick up 
components which oscillate at the Planck scale (see Fig. 
|l|). The precise form of these oscillations depends on 
factor ordering ambiguities (in the coefficients k^ enter- 
ing the constraint) and the use of the Lorentzian (versus 
Euclidean) theory. 

Such a unique wave function generally differs from 
those obtained with boundary proposals of standard 
quantum cosmology. By choosing a real prefactor it is 
always real (for flat spatial slices the evolution equation 
has real coefficients; this no longer holds true for spatially 
curved models) and so cannot coincide with the "tunnel- 
ing" wave function While the "no-boundary" pro- 
posal H also leads to a real wave function, it is imposed 
on the standard Wheeler-DeWitt equation at the Planck 
scale where large deviations to loop quantum cosmology 
occur. Thus, in general its wave function of a universe 
is different from the unique pre-classical solution found 
here. The consistency condition for the initial data in 
loop quantum cosmology may be expressed as sq = 
which is reminiscent of DeWitt's tf)(0) = (!]] (to achieve 
this, an ad hoc Planck potential has been introduced in 
||). However, since these two conditions are imposed 
on completely different evolution equations, the selected 
solutions in general differ. As Fig. |l| shows, there may 
be a good coincidence in certain models, but only if the 
curvature is small at all times which can happen only in 
the absence of matter. 

Contrary to all other proposals for boundary condi- 
tions in quantum cosmology, our dynamical initial condi- 
tions are not chosen to fulfill an a priori intuition about 
the "creation" of a universe but derived from the evolu- 



tion equation which, in turn, is derived from quantum 
geometry, a candidate for a complete theory of quantum 
gravity. Therefore, one equation provides both the dy- 
namical law and initial conditions. As we have seen, the 
critical condition, which crucially depends on quantum 
geometry, emerges from evaluating the evolution equa- 
tion at the state which corresponds to the classical sin- 
gularity. So in contrast to the classical situation where a 
singularity leads to unpredictability, in quantum geome- 
try the regime of the classical singularity fixes ambiguities 
in the wave function of a universe. 

Acknowledgements. The author is grateful to A. 
Ashtekar for discussions and a careful reading of the 
manuscript. This work was supported in part by NSF 
grant PHY00-90091 and the Eberly research funds of 
Penn State. 



* E-mail address: bojowald@gravity.phys.psu.edu 
[1] B. S. DeWitt, Phys. Rev. 160, 1113 (1967). 
[2] C. W. Misner, Phys. Rev. 186, 1319 (1969). 
[3] J. B. Hartle and S. W. Hawking, Phys. Rev. D 28, 2960 
(1983). 

[4] A. Vilenkin, Phys. Rev. D 30, 509 (1984). 
[5] H. D. Conradi and H. D. Zeh, Phys. Lett. A 154, 321 
(1991). 

[6] N. Kontoleon and D. L. Wiltshire, Phys. Rev. D 59, 
063513 (1999). 

[7] A. Ashtekar, Lectures on non-perturbative canonical grav- 
ity (World Scientific, Singapore, 1991). 
[8] C. Rovelli, Living Reviews in Relativity 1, (1998). 
[9] C. Rovelli and L. Smolin, Nucl. Phys. B 442, 593 (1995), 
erratum: Nucl. Phys. B 456, 753 (1995). 

[10] A. Ashtekar and J. Lewandowski, Class. Quantum Grav. 
14, A55 (1997). 

[11] A. Ashtekar and J. Lewandowski, Adv. Theor. Math. 
Phys. 1, 388 (1997). 

[12] M. Bojowald and H. A. Kastrup, Class. Quantum Grav. 
17, 3009 (2000). 

[13] M. Bojowald, Class. Quantum Grav. 17, 1489 (2000). 

[14] M. Bojowald, Class. Quantum Grav. 17, 1509 (2000). 

[15] M. Bojowald, Class. Quantum Grav. 18, 1055 (2001). 

[16] M. Bojowald, Class. Quantum Grav. 18, 1071 (2001). 

[17] M. Bojowal d, Absence of Singularity in Loop Quantum 
Cosmology, gr-qc/0102069| (Phys. Rev. Lett., in press). 

[18] M. Bojowald, The Inverse Scale Factor in Isotropic Quan- 
tum Geometry, in preparation. 

[19] M. Bojowald, Isotropic Loop Quantum Cosmology, in 
preparation. 

[20] T. Thiemann, Class. Quantum Grav. 15, 839 (1998). 
[21] J. M. Sanz-Serna, SIAM J. Sci. Stat. Comput. 6, 923 
(1985). 

[22] A. Ashtekar, J. C. Baez, A. Corichi, and K. Krasnov, 
Phys. Rev. Lett. 80, 904 (1998); A. Ashtekar, J. C. Baez, 
and K. Krasnov, Adv. Theor. Math. Phys. 4, 1 (2001). 



4 



